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Van der Waals explosion of cold Rydberg clusters 
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We report on the direct measurement in real space of the effect of the van der Waals forces 
between individual Rydberg atoms on their external degrees of freedom. Clusters of Rydberg atoms 
with inter-particle distances of around 5 jam are created by hrst generating a small number of seed 
excitations in a magneto-optical trap, followed by off-resonant excitation that leads to a chain of 
facilitated excitation events. After a variable expansion time the Rydberg atoms are field ionized, 
and from the arrival time distributions the size of the Rydberg cluster after expansion is calculated. 
Our experimental results agree well with a numerical simulation of the van der Waals explosion. 

PACS numbers: 34.20.Cf, 32.80.Ee 


The van der Waals (vdW) forces between polarisable 
neutral particles plays an important role in a variety of 
physical phenomena [Ilia. Those forces are weak, ap¬ 
preciable only at very short interparticle distances and, 
therefore, difficult to measure directly, as demonstrated 
by the pioneering experiments on the vdW interactions 
between an atom and a surface ISIE]. Nevertheless, they 
play an important role in weak chemical bindings and 
very often in nanoscale phenomena. The vdW interac¬ 
tion is a key element in the physics of Rydberg atoms 
where, owing to the large distance of the electron from 
the nucleus, the polarizabilities are larger than for atoms 
in the ground state 016 ]. For instance, the Rydberg 
level shifts due to the vdW interactions are responsible 
for the so-called dipole blockade [7] , in which a Rydberg 
excitation is strongly suppressed owing to the presence 
of a previous excitation iHlQj. Experimental determina¬ 
tions of the vdW interaction rely on the measurement of 
observables linked to either internal or external degrees 
of freedom. The first class is represented by the measure¬ 
ments of the level shifts using optical spectroscopy, as in 
an early study m, or in a very recent experiment using 
microwave spectroscopy to monitor the early stage repul¬ 
sion between off-resonantly excited Rydberg atoms m- 
All the above experiments observed the van der Waals 
force indirectly through the effect of the interaction on 
atomic energy levels. The observation of the external de¬ 
grees of freedom of Rydberg atoms was initiated in [13] , 
and progress in that direction is presented in [T4|. Here 
we report the direct observation in real space of the me¬ 
chanical effect of the vdW force on individual Rydberg 
atoms by investigating the formation and subsequent ex¬ 
plosion of Rydberg clusters. Our results enable further 
studies of the vdW interaction as well as the reconstruc¬ 
tion of the spatial arrangement of excited atoms inside a 


Rydberg cluster, in analogy with the experiments on the 
Coulomb explosion of pulse-ionized molecules nans). 

Our experiments, described in caiiH], are performed 
in magneto-optical traps (MOTs) of ^^Rb atoms with 
typical sizes around 150 jam and peak densities of up to 
10^^ cm“^. Two laser beams at 420 nm and 1013 nm 
are used to excite atoms to the 70S Rydberg state, with 
the 420 nm laser detuned by 660 MHz from the 6 P 3/2 
intermediate state (the MOT beams are switched off dur¬ 
ing this excitation). The laser beam at 1013 nm has a 
waist of 110 /im, while the 420 nm beam is focused to 
around 6 /im, thus realizing a quasi-ID geometry in which 
the radial size of the 420 nm laser is comparable to the 
facilitation radius (see below) and hence the excitation 
dynamics takes place predominantly along the x propa¬ 
gation direction of the beam (see Fig. 1). The 1013 nm 
beam lies in the x — y plane at a 45° angle with respect 
to the x-axis, leading to an effective excitation volume 
(defined by the overlap of the two beams) of approxi¬ 
mately 10 “^cm“^. The two-photon Rabi frequency for 
the Rydberg excitation is around 2 MHz. The produc¬ 
tion based on the seed approach (see the following) of 
between one and two Rydberg in the atomic cloud uses a 
short (~ 0.5/is) resonant pulse at small two-photon Rabi 
frequency (0.3 MHz) derived from one of the diffraction 
orders of the acousto-optic modulator controlling the ex¬ 
citation pulses. After excitation and a variable expansion 
time the Rydberg atoms are field ionized and accelerated 
towards a channeltron (see in Fig. 1), where the ion ar¬ 
rival times are recorded and analysed following the ap¬ 
proach introduced in nanni. Our detection efficiency is 
around 50%; for long expansion times the overall signal 
decreases due to spontaneous decay of the Rydberg state. 

In order to be able to determine the size of an expand¬ 
ing cloud of Rydberg atoms over a spatial range of several 
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FIG. 1: Cold atoms in a MOT excited to Rydberg states in an 
elongated interaction volume (continuous ellipse), and subse¬ 
quently expand due to the vdW repulsion. After a variable 
free expansion time they are field ionized by switching on the 
field plates. Finally the ions (blue dots) are electrically accel¬ 
erated through a grid towards a channeltron. The isochronous 
planes defined in the text are schematically represented. 


arrival times are are independent of the ^-position and 
depend linearly on the x- and z- positions, with propor¬ 
tionality coefficients a = —2.31 ±0.08 x 10“^ /is/mm and 
(3 = ±1.34 ± 0.07 X 10“^ /is/mm, respectively. Thus our 
detection process projects the positions of ions created 
in an isochronous plane (which contains the ^-axis and 
makes an angle of 59.9 ± 1.5° with the x — y plane) onto 
a particular arrival time at the channeltron, as shown in 
Fig. 1. We conclude that the ions created in isochronous 
planes separated by a distance As give rise to a temporal 
separation 


At = ± /3^As (1) 

of their arrival times. This dependence was confirmed 
by a measurement of the At delays in the case of MOT 
created varying ^ and keeping x constant and different 
from zero (red circles with black borders in Fig. 2). From 
our data we derive an effective one-dimensional spatial 
resolution of around 9 /im in the single ion regime, where 
the Coulomb repulsion between nearby ions during the 
~ 10 /is acceleration time towards the channeltron can be 
neglected. For the conditions of the experiments reported 
below, this approximation is valid except for the very 
early stages (< 100 /is) of the dynamics, and can be used 
to derive the initial position of the created ions created. 


hundreds of micrometers, we developed a measurement 
technique based on the arrival times of ions at the chan¬ 
neltron. In our apparatus the distribution of arrival times 
of ions originating from the interaction volume is deter¬ 
mined by the geometry of the ion collection apparatus, in 
which the ions are accelerated towards the channeltron 
along a curved path in the x — 2 ;-plane, as schematically 
shown in the upper left corner of Fig. 1. Therefore, the 
ion arrival times depend essentially on their initial po¬ 
sitions along the x and 2 ; axes. The calibration of our 
measurement technique was performed by determining 
the functional dependence of At on the initial positions 
of the ions. To do so, we created small MOTs around 
150 /im in size at different spatial positions using indepen¬ 
dent magnetic fields in the three spatial directions. With 
the MOT lasers switched on, a short (around 1 /is) pulse 
of only the 420 nm beam created on average a single ion, 
in order to avoid aberrations in the collection process due 
to the Coulomb repulsion between the ions. This direct 
ionization was also chosen in order to avoid Rydberg- 
Rydberg interaction effects that could have altered the 
excitation probability and hence the spatial distribution 
of ions inside the MOT. As verified in ED, the direct 
ionization probability is strictly proportional to the local 
atom density. We measured the average ion arrival times 
at the channeltron calculated from 400 experimental runs 
at a given MOT position and thus obtained the depen¬ 
dence of At on the MOT position as shown in Fig. 2. 
To a good approximation, over a range of 0.8 mm the 



FIG. 2: Galibration of the arrival time distributions at the 
channeltron. The center of the arrival time distributions is 
plotted as a function of the MOT position in the x- (blue 
squares), y- (green diamonds) and ^-directions (red circles); 
in each case, the MOT was centred at 0 in the other two di¬ 
rections. The error bars are smaller than the size of the data 
symbols. The dashed lines are linear fits to the experimental 
data within the range (-0.4,0.4) mm. The red circles with 
black borders correspond to a variation in the z-position for 
a fixed x — 460 ± 5 /xm where the arrival times were corrected 
using the x- axis calibration. Their agreement with the x = 0 
data demonstrates that the arrival times are linearly depen¬ 
dent on the X and z positions without cross terms. 
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FIG. 3: (a) Histogram of the probabilities P of ion arrival 
times, Tarr, at A = 80 MHz for two different expansion times: 
130/iS (light blue bars) and 1000/xs (dark blue bars). SB’s 
of P(Tarr) are calculated from Gaussian fits (dashed lines) to 
the histograms, (b) SB, probing the width of the ion cloud 
vs the expansion time for different detunings A: 0 (green 
triangles), 55MHz (blue squares) and 80MHz (red circles). 
On the l.h.s axis SB’s of the arrival times, and on the r.h.s 
axis SB’s of the spatial positions, obtained from the previous 
SB’s by using the inversion of Eq. 0 are shown. The dashed 
lines report numerical simulations of the vdW explosion. 


Two procedures tested in our previous work are ap¬ 
plied to the exploration of the ions created following 
the Rydberg excitation. In order to maximize the van 
der Waals interactions and, therefore, the interatomic 
repulsions, we create Rydberg excitations at small in¬ 
teratomic distances by using the facilitated off-resonant 
excitation explored in [la na na Eg and theoretically 


treated in [23]. In this scheme, following the initial cre¬ 
ation of a Rydberg excitation an additional atom at dis¬ 
tance Vfac from the first one is also excited if the laser 
detuning A matches the Rydberg level shift due to the 
vdW interaction, i.e., hA = Ce/{r^^^) [24|. Facilitation 
implies that close to an already excited Rydberg atom 
all the ground-state atoms located inside a shell of ra¬ 
dius rfac and width Svfac = 1/6 x rfac x S/A (where 6 
is the laser linewidth m) can theoretically be excited. 
The excitation of more than one atom will, however, dis¬ 
tort the facilitation volume due to the additional van 
der Waals interaction between those atoms. For our ex¬ 
perimental parameters, the facilitated excitation of more 
than one atom per facilitation volume is unlikely. For 
the laser detunings of our investigation the values of Vfac 
and Svfac are around 5/xm and 6nm, respectively. The 
first Rydberg atom created through facilitated excitation, 
in turn, facilitates the excitation of another atom, and 
thus a chain reaction is started that continues as long as 
there are ground state atoms inside the next facilitation 
shell. Because the facilitation radius is comparable to 
the transverse dimension of the confining potential the 
cluster of Rydberg atoms created by the chain process is 
extended mainly in one dimension, as shown schemati¬ 
cally in Fig. 4. Since for the large values of A used in 
our experiments the spontaneous off-resonant excitation 
of the initial Rydberg atom needed to start the chain re¬ 
action is small, we use a seed technique based on a single 
resonant pulse leading to the creation of one/two Ryd¬ 
berg excitations within the interaction volume [25l [26] . 
After the creation of the seed a 100 /xs pulse produces 
the facilitated Rydberg excitation. We verify that the 
number of excited atoms reaches a plateau towards the 
end of the excitation time by performing a measurement 
of the mean number of Rydberg atoms by field ioniza¬ 
tion as in [27] immediately after the end of the excita¬ 
tion pulse. From those measurements we derive that for 
A = 55 ± 1 MHz and 80 ± 1 MHz the chain reaction of 
facilitation events produces a quasi-ID Rydberg cluster 
containing 12 excited atoms in both cases. 

We now proceed to use the spatial information con¬ 
tained in the ion arrival times in order to probe the tem¬ 
poral evolution of the size of the Rydberg cloud. After 
the above preparation stage the cluster expands freely for 
up to 1 ms. From our numerical simulations (see below) 
we find that the vdW energy of pairs of Rydberg atoms 
initially created at distance Vfac is fully converted into 
kinetic energy after around 5 /xs, i.e., much less than the 
the 150 /xs lifetime of the 70S Rydberg state. After the 
conversion of the vdW energy into kinetic energy excited 
atoms can, however, decay into lower Rydberg states that 
are not detected by our electric field ionization process, 
which only field ionizes and detects Rydberg atoms with 
principal quantum number n > 50. By repeating the ex¬ 
citation/detection cycle 400 times we determine the spa¬ 
tial distribution of the Rydberg cloud. Here we assume 
























4 



x(|jm) 

FIG. 4: Simulated van der Waals explosion in the x — z plane 
of a ID cluster of Rydberg atoms (white spheres represent 
the initial positions of the Rydberg atoms, and black spheres 
indicate the positions after 60 /is). Each segment (black and 
white) of the trajectories corresponds to 5 jis. 


that the spontaneous decay of the Rydberg excitations 
occurs independently for each atom, so that the arrival 
time distribution of the non-decayed subset of atoms re¬ 
flects that of all the initially excited atoms. Typical his¬ 
tograms of the arrival time distribution are reported in 
Fig. 3(a). We measure the standard deviations (SB’s) of 
those distributions as function of the free expansion time, 
as shown in Fig. 3(b) (left scale). By inverting Eq.Q, 
from the arrival time SB’s we derive the spatial width 
SB’s as a function of that time, as shown in Fig. 3(b). 
For both values of A different from zero the vdW explo¬ 
sion produces a large spread in the arrival times. In a 
test experiment at A = 0 (with 12 excitations on aver¬ 
age), for which no vdW explosion should occur, the cloud 
expansion at 1 ms is small and compatible with the ex¬ 
pected expansion due the thermal motion. As expected 
from the strong r-dependence, the van der Waals force 
acts mainly at short interatomic distances and, therefore, 
the expansion of the cloud is essentially ballistic over the 
full temporal range. In fact, the large initial acceleration, 
^ 5 X 10^ ms“^ at A = 80 MHz and Vfac^ reduces to one 
percent at 2rfac‘ At that detuning the theoretically ex¬ 
pected overall increase in velocity for two Rydberg atoms 
with a corresponding initial vdW energy is « 600 jxm/ms 
compared to the average thermal velocity of MOT atoms 
of ~ 120 /rm/ms. 

Fig. 3 (b) also shows the results of a numerical sim¬ 
ulation of our experiment (see Fig. 4), which includes 
the van der Waals repulsion, the thermal motion of the 


atoms and the Coulomb repulsion after field ionization, 
which acts during the ^ 10 /is time-of-flight towards the 
channeltron. The creation of the clusters is simulated by 
placing Rydberg atoms at a distance Vfac with a probabil¬ 
ity weighted by the square of the Rabi frequency (since 
for our experimental parameters the excitation process 
is incoherent [28]) that varies perpendicularly to the x- 
direction. The constraint in the interparticle distance 
together with the spatially varying excitation probabil¬ 
ity gives rise to quasi-ID clusters that grow close to the 
x-axis with small (of the order of the width of the 420 nm 
beam) and randomly distributed radial separations from 
it. In our simulations the nuclear motion is treated clas¬ 
sically with an adiabatic evolution of the electronic states 
and neglecting all nonadiabatic couplings [29]. For long 
expansion times of the clusters. Fig.3 (b) shows excel¬ 
lent agreement between the experiment and the numeri¬ 
cal simulation without any adjustable parameters. Below 
200 fis the simulation slightly overestimates the effect of 
the Coulomb explosion, which may be due to our incom¬ 
plete knowledge of the exact trajectories of the individ¬ 
ual ions. In the simulation a linear acceleration due to a 
perfectly homogeneous electric field along the x-axis was 
assumed, which is probably not an accurate representa¬ 
tion given the geometry shown in Fig. 1 (these deviations 
will be investigated in detail in future work). The simu¬ 
lation also reveals that the deviation of the clusters from 
a perfect ID configuration leads to a smaller expansion 
(by about a factor of 2) in the x-direction than in the 
radial directions. For hypothetical radial beam widths 
smaller than around 2 /im, by contrast, the simulations 
show a larger expansion along the x-axis. 

In summary, we have measured the mechanical effect 
in real space of the vdW interaction between Rydberg 
atoms excited at close range. These measurements were 
made possible by interpreting the arrival time distribu¬ 
tions of ions at the channeltron after field ionization of 
the Rydberg clusters. Our observations indicate that 
the facilitated creation of Rydberg excitations at large 
values of the detuning, which leads to small separations 
between the excited atoms and hence to a large inter¬ 
action energy, makes it possible to observe the van der 
Waals repulsion even in magneto-optical traps, without 
the need to go to lower temperatures such as those com¬ 
monly reached in Bose-Einstein condensation. Whilst the 
present experiment, which projects the 3D Rydberg cloud 
onto a single time dimension, cannot observe details of 
the vdW interaction, that resolution can be achieved in 
future experiments with 2D or 3D resolution based on 
large area imaging, precise multi-hit timing information 
and high detection rates, such as in the helium BEG ex¬ 
periment [30]. Also, the Rydberg atoms accelerated as 
shown above could be used for controlled collision ex¬ 
periments. Finally, our results demonstrate the intrinsic 
mechanical instability of off-resonantly excited Rydberg 
gases, which has implications for proposed applications 
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in quantum computation and quantum simulation m- 
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